Two-point correlation function of the fractional Ornstein-Uhlenbeck process 
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We calculate the two-point correlation function {x{t2)x{ti)) for a subdiffusive continuous time 
random walk in a parabolic potential, generalizing well-known results for the single-time statistics to 
two times. A closed analytical expression is found for initial equilibrium, revealing non-stationarity 
and a clear deviation from a Mittag-LefHer decay. Our result thus provides a new criteria to assess 
whether a given stochastic process can be identified as a continuous time random walk. 
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I. INTRODUCTION 



During the last decades, continuous time random walks (CTRWs) have been used to model a wide variety of 
stochastic processes which can not be described by the usual type of Langevin equations. Especially, CTRW statistics 
have been successfully applied to systems exhibiting non-Gaussian probability distributions. For a review we refer 
the reader to the comprehensive review articles by Metzler and Klafter [H [2]. It has become evident however, that 
due to their non-Markovian character, assessing the probability distributions at a single time is insufficient in order 
to determine the nature of the underlying stochastic process. To this end, multiple time probability distributions 
of CTRWs have been investigated in a recent line of research [2J 0] |5J IHJ IZJ |S] . These probability distributions can 
be used to determine multiple time correlation functions, which in turn allow for a comparison with experimental 
data. Experiments in which such correlation functions play a crucial role have recently been reported in the context 
of protein conformational dynamics [HI HO]) and blinking nanocrystals [IT]. 

The purpose of the present paper is to provide the two-time correlation function for a subdiffusive CTRW-process 
of a particle in a harmonic potential, also referred to as fractional Ornstein-Uhlenbeck process. A general form of 
this correlation function as well as a closed analytical expression for initial equilibrium position of the particle are 
derived as our main results. It is clear from many textbook examples in physics, that motion in a parabolic potential 
is of utmost importance, and has been used as an approximation e.g. in [S^. Our results could therefore be readily 
compared with data from existing and future experiments. 

The remainder of this paper is organized as follows. In the next section we present the basic equations for the 
CTRW model. Then, well-known results for the single-time statistics are rederived. The two main sections follow, in 
which we derive our two main results, the general expression for the correlation function and a special analytic form. 
We conclude with a few remarks. 



II. COUPLED LANGEVIN EQUATIONS OF THE FRACTIONAL ORNSTEIN-UHLENBECK PROCESS 

In this section we briefly review a well-known representation of CTRWs in terms of coupled Langevin equations 
[T2j. In the continuum limit of infinitesimal step lengths, the one dimensional motion of a random walk particle in a 
parabolic potential V{x) = 1/2 'fx^ is described as: 



dX{s) 
ds 

dt{s) 
ds 



-jX{s) + r,is), (1) 
r{s). (2) 



In this framework the CTRW is parametrized by the continuous path variable s, which may be regarded as arclength 
along the trajectory. X denotes the physical space and t the physical time. Both are given as stochastic processes 
in the 'eigentime' s. Their statistics are determined by the properties of the stochastic variables 77(5) and t(s). We 
only consider the special case of statistically independent increments rj and r, corresponding to uncoupled jump 
lengths and waiting times. Assuming 77(5) as a standard Langevin force with properties {ri{s)) = and (fy(s)?7(s')) — 
2aS{s — s'), we see that X(s) accordingly corresponds to the ordinary Markovian Ornstein-Uhlenbeck process |13j . 
CTRW characteristics enter via the process t{s). Its increments r(s) are assumed to be broadly distributed such 
that t{s) represents an asymmetric Levy-stable process of order a with < a < 1. Levy-stable processes of this 
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kind induce a diverging characteristic waiting time (Ms)). As a result, the second moment of the corresponding free 
diffusion process, i.e. setting the force term in Eq. to zero, yields subdiffusive scaling (a;^(i)) ^ t". Therefore we 
refer to the stochastic process specified by Eqs. ([T]) — ^ as a subdiffusive CTRW in a parabohc potential, or fractional 
Ornstein-Uhlenbeck process [Tj . 

The CTRW in physical time t is determined by the subordinated process X{t) = X{s{t)), where s{t) is an inverse 
Levy-stable process jT4]. We define the probability distributions (pdfs) of X{t) as 

f{x2M\XlM) = {5{X2- X{S2))5{S2~ s{t2))5{x^- X{si))5{s^- s{h))), (3) 

and accordingly for multiple times. It is a basic result, that these n-point pdfs can be expressed as an integral 
transformation [5] 

/■oo /"OO 

f(.{xi,U}) = / dsi... dsn h{{si,U})fM{{xi,Si}), 
Jo Jo 

(4) 

here h({si,ti}) denotes the n-point pdf of the process s(t) and fMii^iT Si}) the n-point pdf of the Markovian process 
X{s). Equation Q states that the pdfs of the fractional Ornstein-Uhlenbeck process are determined by a trans- 
formation of the corresponding Markovian pdfs. The integral kernel is generated by an inverse Levy-stable process. 
From this integral transformation it is obvious, that all moments of the process X{t) are obtained by averaging the 
moments of the Markovian process X{s) with the subordinator s{t). In the single and two point case, h assumes a 
simple analytical form in Laplace space. 



III. SINGLE-TIME MOMENTS 

Let us first consider the well-known single-time moments {x{t)) and {x^{t)). In this case Eq. Q leads to 

/•oo 

/ ds{x"-{s))h{s,t). (5) 
Jo 

Here the single-time pdf h{s,t) is given as [TS] 

(6) 

La denotes a single-sided Levy distribution of order a. As mentioned above, calculations involving h are more 
conveniently performed in Laplace space. We define Laplace transforms in the usual way: C{g{t)} — g{t)e~^^dt, 
accordingly for multiple variables tj A^. Functions with Laplace space argument always denote the Laplace 
transform if not otherwise indicated: g{\) = C{g{t)}. In Laplace space h{s,t) reads h{s,X) = A"~^e~^ * and with 
{x{s)) = xoe"''"* (see e.g. Ch. 3 in P2|) we can write the Laplace transform of {x{t)) 

poo A"^i 
{x{X)) ^ xo d.A"-ie-(^+^")^ = xo— — . (7) 

Jo "f -t A 

Performing a series expansion of 1/(7 + A") and term by term Laplace-inversion, the final result can be expressed 
with the help of the one-parameter Mittag-Leffler function Ea [TB] 

^-'{^}-tP^:-^^'-yn. (s) 

I 7 + A" J ^r(na+l) 

V y n— 

and therefore 

{x{t)) =xoEa{-7n. (9) 

Likewise we can calculate the second moment {x'^{t)). Note that (a;^(s)) = {xq — cr/7) exp(— 27s) + cr/j ([13]) Ch. 3), 
and we obtain the result 

{x\t)) = (^xl - Ea{~2^n + ^. (10) 

In the long-time limit the mean-squared displacement approaches the thermal equilibrium value {x'^)eq ^ <^ll = 
k^T 1^. Above results for the single-time moments have been derived in the context of the fractional Fokker-Planck 
equation [T7] . 
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IV. CORRELATION FUNCTION 



The calculation of two-time moments is more involved but follows the same procedure. In the following we restrict 
our consideration to the simplest case {x{t2)x(ti)), higher order moments can be obtained along the same lines. With 
the help of the integral transformation Eq. d4|, we can express the correlation function as before 



ds2 / dsi (x{s2)x{si))h{s2,t2:,si,ti). 



(11) 



Here, the correlation function {x{s2)x{si)) of the standard Ornstein-Uhlenbeck process is well known [13^ and follows 
from the Langevin equation ([T]) in a straightforward way 



{x{s2)x{si)) = fa;^_-V"^(^2+«i) + ^g-7N2-sil)^ 
V 7/ 7 



(12) 



On the other hand, the two-time pdf h{s2,t2\ si,ii) has a closed analytical form in Laplace space only [B] 



ft-(s2, A2; Sl, Ai) = 



1 



92 



0(s2 - si)e 



-(Ai+A2)°Sl-A2 (S2-Sl) 



A2A1 ds2dsi 



(13) 



Performing the integrations in Eq. (11 1 with these two expressions, we obtain the following result for the correlation 
function in Laplace space 



(x(A2)a;(Ai)) = a(Ai -I- A2)-" 



1 A° 



1 A? 



Ai(7+AJ) A2(7 + A?) 



° 7/ AiA2(27+(Ai + A2)") 



7 + A^ 7 + A? 



+4 



1 x'^-^ 



1 Ar^ 



1 



_A2(7+A?) Ai(7 + A^) AiA2_ 
Laplace inversion yields our first main result, the correlation function in physical time 



(14) 



[E^{-jt^) + i?«(-7i?)] +[xl--]Zsh,t2n,ti) 



(.(t2).(t0) = -[9^2+91; 

+ -[Eo.i-jt^) + E^i-jt?) - 1]. (15) 
7 

Here we have introduced the two-time fractional integral operator {dt^ + c?ti)~", which can be represented as [B] 



_d_ 

df2^dh 
and whose Laplace transform reads 



g{t2,h) 



r(«) Jo 



Min{t2,tx) 



dtt°'^^g(t2 -t,ti- t), 



(16) 



d_ d_ 



g{t2,h) ) = (A2 + Ai)-"5(A2,Ai). 



Furthermore, Eq. (15 1 contains the Laplace transform of the two-time pdf h of the process s{t) 

= ds2 dsie^"^^^-"i^i/i(s2,i2;si,ti), 
Jo Jo 



(17) 



(18) 



evaluated at U2 = ui = 7. The occurence of this Laplace transform is simply a consequence of the first term in 



Eq. (12l and the transformation Eq. (11). In the following we use the properties of the pdf h(s2,t2', Si,ti) in order to 
validate certain limit cases of Eq. ( 15 
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A. The special case ti = 



We first rederive the autocorrelation function {x{t)x{0)) from Eq. (15). For ii = we know that /i(s2, t2; si, 0) ~ 
h{s2,t2)S(si) and thus Zs(7,t2;7,0) = .^3(7,^2) = ^'c((— 7*2)- The fractional integral operator has zero contribution 
for ti — 0. We therefore recover the well known result (see Eq. (|9l) 



{x{t)x{0)) ^ xlE^{~-ie) = XQ{x{t)). (19) 



B. The limit t2 ti 



For ^2 ~^ ti Eq- (15 1 should reduce to the quadratic moment Eq. (10 1. In this limit the two-time fractional integral 
reduces to the familiar Riemann-Liouville fractional integral 

d d 



tl^i \dt2 ' dti 



Eo.{'lt1) = oA";"^a(-7t?) 



1 



(20) 



r(«) Jo 

Using the series representation of the Mittag-Lefflcr function, the convolution integral of each term is calculated in 
Laplace space in a straightforward way 



00 



7 ^ Ai 

n— 1 



(21) 



Therefore 



oi^i-"i?„(-7i?)= -(l-i?a(-7i?)) 



(22) 



Further note that limt^^ti /i(s2, ii; Si, ^i) = 5{s2 — Si)h{si,ti). A quick calculation then yields ^3(7, ti; 7, ii) — 
Ea{—2jtf). Substituting this expression and Eq. (22l into Eq. (15l leads to the result limt^^tj (a:(t2)a:(ii)) = {x^{ti)). 



The limit a 



In the limit a — > 1 the first moment of the CTRW waiting time distribution is finite and we expect that Eq. (15 1 
yields the standard correlation Eq.(12). This can be seen as follows 

t2)S{si 



a — » 1 limit: lim^^i h{s2, ^2! Si,ti) — 6{s2 — ^2)^(^1 — ii) (see Eq. ( 13 1), and thus lim^^i ^^(7, i2; 7j ^i) = e 
In the same limit the Mittag-Leffler function reduces to the usual exponential function and the two-time fractional 



The pdf h factorizes into delta-functions in the 

7(t2+tl) 



integral to a normal integral. The result of this calculation is 



_d_ 
dt2 



d_ 



-7*2 



e-^*i) = Q{t2-h)- (1 



-7(t2-tl) 



+Q{h-t2)- fl + e-^(*i-*^) 



-7*2 



= -7*2 



-7*1 



= -7*1 



(23) 



Substituting these expressions into Eq. ( 15 1 yields the standard correlation of the Ornstein-Uhlenbeck process Eq. ( 12 1 
as expected. 



V. INITIAL EQUILIBRIUM POSITION 



If we choose as special initial value the thermal equilibrium displacement {x'^)eq, we can determine a closed analytical 
expression for the correlation function {x{t2)x{ti)). For Xq ~ (7/7 Eq. (151 leads to 



(x{t2)x{ti))Q 



d_ 

dr2 



d_ 



[E^(-jt^) + E^i-jf^)] 



+ -[E^{—ft^) + E^{-^t'})-l]. 
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(24) 



Time tj 



FIG. 1: The correlation function Eq. (28l for j = a — 1, t2 > ti = 1 and three different Q-values: a = 1.0 (dashed line), 



a = 0.8 (solid), and a = 0.5 (dotted). For a — 1 the decay is exponential, see Eq. ([29 




Time U 



FIG. 2: The correlation function Eq. (28 1 for j — a — 1, t2 > ti = 1 and a = 0.8 compared with the Mittag-LefBer function 



Ec,{—{t2 — ii)") for same parameter values. 



Here, the two-time fractional integral reads explicitly 



d_ d_ 



[^a(-7i2)+£^a(-7i?)] = e(t2-tl) 

+e(ii -ia) 



r(a) Jo 



1 



dtr-^E^{-^{t2-tr) 



r(«) Jo 

The second integral can be evaluated with the series representation of the Mittag-Leffler function. We note that 

/ dii"-i(i2-i)""=^^2i^i a,-an,a + l;-i , (26) 

Jo "V ^2/ 

where 2^i(a, ^, c; 2;) denotes the Gaussian hypergeometric function [18^. Consequently 



r(a) Jo 



dtt"-'E^i-jih-tr) 



ti (-7^2)" ^ . 

TT773^^2i^i ( a, -ari, a + 1 



r(a + l) ^ r(an+ 1) 

^ n— 



(27) 
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Upon substituting this expression as well as Eq. ( 22 1 back into Eq. ( 25 1 , we obtain our second main result from 
Eq. (24), the correlation function for initial equilibrium in closed form {t2 > ii) 

{x{t2)x{ti))o 



r(a + 1) ^ r(Q;n + 1) V t2/ 7 



The limit cases hold as before, as one can easily check. In particular for a — > 1 {t2 > ti) 



lim{x{t2)x{ti))o = -e' 

a— >1 ^ 



-7(t2-tl) 



(28) 



(29) 



With above equations one can see that the CTRW correlation function becomes a function of the time difference 
t2 — ti only in the limit a —^ 1. This characteristic property is not evident from the autocorrelation function Eq. ( 19 1. 



Eq. ( |28| therefore clearly demonstrates, that the fractional Ornstein-Uhlenbeck process is a non-stationary process 
for diverging characteristic waiting times (i(s)) of the random walk particle. 

Fig. [l] shows the correlation function Eq. ( pS]) for different a-values. In Fig. [2] the deviation from a Mittag-Leffler 
decay is clearly visible. The two curves in Fig. [21 only agree for ti = 0. 



VI. CONCLUSION 



We have derived expressions for the correlation function of the fractional Ornstein-Uhlenbeck process, starting 
from a representation in terms of coupled Langevin equations. For initial equilibrium position of the random walk 
particle we were able to obtain the correlation function in analytical form. This correlation function exhibits two 
main features. Firstly, it demonstrates the non-stationarity of the CTRW in a harmonic potential for < a < 1. 
Secondly, it clearly deviates from the familiar Mittag-Leffler decay. Both findings are in contrast to the properties of 
the autocorrelation function Eq. ( 19 ), to which our result converges for ti — 0. Considering the widespread application 
of the single-time quantity Eq. (19), comparing experimental data with the correlation function Eq. (28) should be a 
worthwhile task. In the end only multi-time statistics are appropriate measures in order to distinguish CTRWs from 
other non-Mar kovian processes. 
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